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Abstract 

Near-soliton scanning light-beam propagation in media with both 
delayed-response Kerr-type and thermal nonlinearities is analyzed. 
The delayed-response part of the Kerr nonlinearity is shown to be 
competitive as compared to the thermal nonlinearity, and relevant 
contributions to a distortion of the soliton form and phase can be 
mutually compensated. This quasi-soliton beam propagation regime 
keeps properties of the inclined self-trapped channel. 

1. Introduction 

A competition between the cubic and thermal medium nonlinearities un- 
der different laser pulse durations seems to be firstly considered in [IJ. Re- 
cently we noted some special features of the soliton planar (two-dimensional) 
light beam propagation in a nonlinear medium with both the promt Kerr- 
type nonlinearity and the inherently delayed-response thermal one. Under 
condition of scanning a powerful light beam, the Kerr-type nonlinearity was 
shown to be able to become dominant, while the thermal one being a per- 
turbation exerting an influence upon the form and phase of a Kerr soliton. 
In this situation the planar beam keeps the auto-waveguide properties. 

In this communication, we pay attention to some effects which occur when 
taking account of the material response of both simultaneously existing non- 
linearities. Fast scanning light beam is the case where the dynamic cubic 
nonlinearity can be appreciable. It is shown that under certain condition the 



delayed-response part of the Kerr-type nonlinearity and the thermal nonlin- 
earity can be mutually cancelled. 



2 



2. Formulation of the problem 

A system of equations governing the planar light beam propagation, al- 
lowing for both above nonlinearities and under the assumption of slight ab- 
sorption, reads ||. [4|] 



E Z + -E t + ^E xx = -i-(An^ + An^)E, An™ = n T T, 
u 2k no 

r(An^) t + AnW = n 2 \E\\ Z = X T XX + ^\E\ 2 . (1) 

27TpCp 

Here E is the beam field, T is the medium temperature measured with respect 
to the equilibrium one T , 7 is the linear absorption coefficient, x is the 
temperature conductivity coefficient (the heat conductivity coefficient o = 
c pPX)i c p is the specific heat capacity at constant pressure, p is the medium 
density, no is the linear refraction coefficient, n 2 is the Kerr constant, n = 
n Q + An, An = An^ + An^ characterizes the contributions of the Kerr- 
type nonlinearity (K) and the thermal one (T), respectively, ko is the wave 
number, k = n k Q , r is the relaxation time of the Kerr-type nonlinearity, X 
is the lateral coordinate directed along the interface, Z is the longitudinal 
coordinate directed into the bulk of the nonlinear medium, u is the group 
velocity in the first approximation of the dispersion theory, t is the running 
time. There is no Y dependence according to the two-dimensional model. 

Let the action of a powerful radiation on the nonlinear medium be realized 
via translational scanning along the X axis of the planar beam, the incidence 
direction being normal to the interface and coincident with the Z axis. Also, 
let v be the scanning velocity and D be a width of the beam. Assume that 
the field E and the temperature T are dependent only on the moving frame 
coordinates X' = X — v(t — (Z/u)) and Z' = Z, being independent on 
time explicitly, i.e., the stationarity condition holds in the moving reference 
system. In this system we have 

Ez , + \E X , X , = -i—(AnW + An^E, An^ = n T T, 
2k n 

An« « n 2 \E\ 2 + rvn 2 (\E\ 2 ) x ,, -vT x , = ^^\E\ 2 . (2) 

OTlpCp 
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The expression for An^ is obtained from an approximate solution of the 
corresponding equation (1) under the assumption of smallness for r ||. Note 
that the small term containing the derivative of \E\ 2 can have another origin 
related to the Raman process |J. Besides, the second order derivative in 
temperature conductivity equation (1) may be neglected when (x/vD) « 1 
[0], which permits to represent T as T ~ \E\ 2 dX' . Eventually, using (2), 
we obtain a resulting equation for E representable in a convenient dimen- 
sionless form: 

%A Z - A xx - 2\A\ 2 A = a(\A\ 2 ) x A ± (3A T \A\ 2 dx . (3) 

J — oo 

Here the following dimensionless quantities are introduced: 



D 2kD 21 4 \7iaTf 



( D\ n 2 QO 3 (rT Q Dn 2 T jj.8 



g = 16irk 6 crT — — , a = 32wk 6 — 9 v, (3 = —, (4) 
\n J c cn z g g 

T 

5 = 2kDxiv~ 1 , n = 2(kD) 2 — \n T \. 

n 

The signs ± for the last term in (3) correspond to the defocusing or focusing 
effect of the thermal nonlinearity, respectively. 

The coefficients a and f3 are seen from the derivation of the system (1) to 
be small (a << 1 and (3 « 1), which allows Eq. (3) to be considered as a 
perturbed nonlinear Schrodinger equation (NSE). In the following, we shall 
conceive that the soliton relevant to the non-perturbative NSE is formed at 
the interface. Our task is to describe a behavior of the soliton under the 
action of the perturbation (3). For definiteness, we restrict ourselves to a 
consideration of the defocusing thermal nonlinearity. 



3. Results. 

In the adiabatic approximation where a distortion of the soliton form is 
neglected, we seek a solution of Eq. (3) as 

A s (x, z) = 2r\{z) exp(—i9(z))sechy, 
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where y = 2r\{x — k(z)), 9 = [£(z) /i](z)]y + A(z) and z-evolution of the 
parameters £, 77, k and A is given by the formulae 0: 



£ z = -3? / .Rtahm/ sechj/ e* 9 cfa/, r^ z = — —Of / Rsechy e ie dy, 

Z J — oo Z J — oo 

k z -^=—^s J^yR sechy (5) 



A z - 2£k z + 4(£ 2 - r? 2 ) = —3? / #(1-2/ tahny)sechy e i9 dy. 

2n J-oo 



2rj 

Here the perturbation is given by 

R = a{\A\ 2 ) x A + (3A ( X \A\ 2 dx'. (6) 

J — oo 

Solving the equations (5) we find 

£ = £o + jzJVoZ, V = Vo, « = 4 £o z + -^7 rj 2 z 2 + « , 
15 15 

A = 4(£ 2 + v l - lp v )z + Vl z 2 + |(^7 % 2 ) 2 + Ao, (7) 

where 

7 = 16a rj 2 - 5/3. (8) 

The soliton amplitude 2rj is seen to be independent on the perturbation while 
an angle of inclination of the self-trapped channel to the z axis (related to 
the parameter £) varies, an increase or decrease of this angle being deter- 
mined by the sign of a certain combination 7 (8) of the parameters a and /3. 
The function k(z) characterizes a soliton centre location and depends on the 
perturbation also by means of 7. The soliton phase 8(z) acquires a nonlinear 
dependence on z, the latter being completely determined by perturbations 
nonlocal in x. 

Finding corrections to the soliton form enables a conclusion to be drawn in 
relation to keeping or non-keeping a soliton localization under a perturbation. 
In the framework of the first-order approximation where a solution of Eq. (3) 
is written as A — A s +Ai, the correction Ax is evaluated by the formula [0. [| 

1 ( r ^ w 
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+r]' z e _2ie sech^ 



dX- 



6(A) 



2i\x 



-oo (A-A 1 )(A-A 1 ) /' 

where Ai = £ + ir], the bar implies complex conjugation, and 6(A) depends 
on the perturbation (6) (see p|): 



ITX 



b(z, A) = -(A-£) 



12a 
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(A_A 1 )(A-A 1 ) 



e *(A-2A K ) sech 



7rA-£ 

2 r? 



Restricting ourselves to finding asymptotics, we obtain 



x 



+oo : 



x — > — oo : 



y-W 



16 2 1 / 

— «?/ + — 72/(2/ - 2) 

15 15 

16 2 1 / 



Thus, the soliton localization keeps in the presence of the perturbation (6) 
though a slight disturbance of the initial soliton x-symmetry is observed. 



4. Discussion 

To illustrate the plausibility of the idea proposed, let us adduce some 
typical parameters of available nonlinear condenced media and optical beams 
(order-of-magnitude estimates): 

n ~ 1, n 2 ~ 1(T 13 . . . 1(T 12 CGSE, 

r ~ 10" 12 ...10~ n s, L)~10 _1 cm, v ~ 10 8 cm/s, 

pc p ~ 10 7 erg/cm 3 K, 7 ~ 10~ 4 cm _1 , k ~ 10 5 cm~ 1 . 

The magnitudes of the coefficients a and f3 calculated with the above pa- 
rameters confirm our conjecture concerning the smallness of a and (3. Thus, 
the pertinent experimental conditions seem to be quite realizable, including 
a corresponding scanning velocity value. The last is achievable by means of 
up-to-date deflection techniques. 

The perturbation (6) is described by two terms each of them reflects 
nonlocal character of the delayed-response nonlinearity. The non-locality of 
the perturbation determined by the inertial Kerr nonlinearity is of differential 



dn 



dT 



10" 4 ...10" 3 K- 1 , 



6 



type while the non-locality due to the thermal nonlinearity is associated with 
the integral type. In the moving reference frame the above non-locality is 
actually similar to the spatial dispersion influencing the soliton stability [|]. 

If the coefficients a and (3 are of different signs, the action of both per- 
turbations leads to a summation effect. On the other hand, if a and (3 are of 
the same sign, the nonlinear perturbations able to be compensated partially 
or completely. This is seen from relevant corrections (7). A complete com- 
pensation occurs at 7 = 0. Let us write down a and j3 isolating the explicit 
dependence on the scanning velocity: a = a*v and (3 = (3*v _1 (expressions 
for a* and f3* can be readily obtained from (4)). Then we get from the 
condition 7 = 0: 

« 2 = («) 



Hence, a feasibility of nonlinear distortion compensation is due to differ- 
ence in v dependence of both nonlinearities. Namely, the influence of the 
Kerr-type nonlinearity increases with increasing v, while that of the thermal 
nonlinearity decreases. Note that the relation (9) resembles the resonance 
condition for the oscillating contour. 

Thus, a compensation of distortions caused by the delayed-response part 
of the cubic nonlinearity and defocusing thermal nonlinearities leads, in the 
adiabatic approximation, to the "pure" soliton having the envelope of sech- 
type and the phase linear in z. In this case, there remains a linear correction 
due to the thermal nonlinearity. It is worth to point out that one can also 
use tilted pulses M which resemble translationally scanned beams. 



* * * 
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